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Abstract: The aim of this paper is to study the law of the last passage time of a linear diffusion to a 
curved boundary. We start by giving a general expression for the density of such a random variable under 
some regularity assumptions. Following Robbins & Siegmund, we then show that this expression may be 
computed for some implicit boundaries via a martingale method. Finally, we discuss some links between 
first hitting times and last passage times via time inversion, and present an integral equation (which we 
solve in some particular cases) satisfied by the density of the last passage time. Many examples are given 
in the Brownian and Bessel frameworks. 
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1 Introduction 

1.1 Motivation 

Let I G [— oo,+oo[ and consider a linear regular conservative diffusion (X t ,t > 0) taking values in 
/ = (£, +oo[ with +oo a natural boundary. 

Let / : [0,+oo[ — > [I, +oo[ be a continuous function and define C(/) : = inf{t > 0; f(t) = £} G]0,+oo]. 
The aim of this paper is to compute the law of the last passage time of X to the boundary / before time 
C(/): 

G f := su P {0 < t < C{f) ; X t = f(t)}, 

(=0 if{o<t<c(/);*t = /(*)} = 0)- 

This problem was essentially addressed (in a far more general framework) in the literature through the 
study of additive functionals. We refer in particular to Getoor & Sharpe [GS73, Proposition 3.3] where 
the law of the last exit time from a Borel set D of a general Markov process is computed, with the help of 
the potential kernel of an additive functional associated to D. In [Por67], the case of R"-valued symmetric 
stable Levy processes is tackled, and the law of the last exit time from a Borel set D C W 1 is obtained in 
term of the equilibrium measure of D (see also Takeuchi [Tak67] for a study of moments). 

In the set-up of diffusions (which is our concern), one of the main result is due to Pitman & Yor 
[PY81], who compute the law of Gf when f(t) = a is a constant boundary and (X t ,t > 0) is a transient 
diffusion going to +oo a.s. (Another proof is given in [PRY10, Chapter 2, p. 38] using a formula for the 
last passage time of a continuous local martingale to a constant boundary.) We shall recover their result 
in Example 2 below. 

On the other hand, much emphasis has been put, quite naturally, on the study of the first passage 
time of X to the boundary /: 

T f :=inf{t>0;X t = /(*)}, 

(=+oo ii{t>O;X t = f(t)} = 0). 
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Many ideas have emerged to solve this problem: changes of measure and partial differential equations 
[Sal88, Gro89], integral equations [RSS84, BNR87, GNRS89], martingale methods [RS70, Nov81], inver- 
sion of time [PY81, AP10]... We shall see that some of these methods may be applied to get information 
on last passage times. 



1.2 Notations 

Let (X t ,t > 0) be a linear regular diffusion taking values in / = (£, +00 [ with +00 a natural boundary. 
We assume that (X t ,t > 0) is conservative, in the sense that it has an infinite life-time, but we make no 
assumption on the nature of the boundary point £. Let V x and denote, respectively, the probability 
measure and the expectation associated with X when started from x > I. We assume that X is defined 
on the canonical space il := C(R + — s- I) and we denote by (J-*, t > 0) its natural filtration. 

Wc denote by m(dx) — p(x)dx its speed measure, which is assumed to be absolutely continuous with 
respect to the Lebesgue measure and by s its scale function, which we assume to be of C 2 class. With 
these notations, the infinitesimal generator of (X t ,t > 0) reads: 

d 2 

G= -, for (t, x) G R + x I. 

p(x)ox s'(x)ox 

It is known from Ito & McKean [IM74, p. 149] that (X t ,t > 0) admits a transition density q(t,x,y) with 
respect to m, which is jointly continuous and symmetric in x and y, that is: q(t,x,y) = q(t,y,x). In 
the remainder of the paper, we shall always assume that q is a C 1 ' 2 ' 2 class function on ]0,+oo[x/ x /. 
In particular, for any y £]£, +oo[, the function (t,x) 1 — > q(t,x,y) is a solution of the partial differential 
equation: 

5q = ^ on]0,+oo[x]£,+oo[. (1) 

1.3 Organization of the paper 

• We start in Section 2 by giving, under some regularity assumptions, a general expression for the density 
of the r.v. Gf. 

F x (G f G dt) = $(t)q(t, x, f(t))dt, (0<t< C(/)) (2) 

where 

Observe that the dependence on the initial state only appears through the transition density. Several 
examples of application of this formula are given, involving Brownian motion and Bessel processes with 
linear, squared root and square boundaries. 



• In Section 3, we follow Robbins & Siegmund [RS70] and present a martingale method for computing 
the function $ which appears in the previous formula. This gives us the law of Gf for a large class of 
implicit boundaries /. 

• We then discuss, in Section 4, some relations between first hitting times and last passage times through 
time inversion. 



• Finally, for all initial states x such that ¥ x (Gf > 0) = 1 (i.e. the law of Gf under P x has no atoms), 
Formula (2) implies that: 

C(/) 

f>(t)q(t,x,f(t))dt = l. 
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We shall see in Section 5 that in some cases, this relation characterizes uniquely the function hence 
the law of G/ . 



2 The density of G f 

We start by establishing a general formula for the density of the last passage time Gf. Let (0 t ,t > 0) 
denote the translation operator defined by: 

Vu>o, (foe t )(u) = f{t + u). 

In the following, we shall distinguish two cases, depending on whether the process (X t ,t > 0) remains 
above the boundary / after the last passage time Gf (we shall say that / is a lower boundary) or under 
the boundary / (resp. / is an upper boundary). 

2.1 Lower boundaries 

Let / : [0,+oo[ — > [£, +oo[ be a continuous function, which is of C 1 class on ]0,+oo[ and define £(/) := 
inf{i > ; f(t) = £} G]0, +oo]. If £(/) = +oo, we suppose (to ensure that Gf ^ +oo a.s.) that: 



Va; 6 /, 



lim X t - f{t) > = 1. 
*-K(/) 



(3) 



This implies in particular that, for any value of C(/) G]0, +oo], the diffusion (X t , t > 0) remains above the 
boundary / on the time interval ]G/, C(/)[- 




C(/) = +^ 




Gf C(f) ^ 
Let H be the function defined by: 

H ■ (t, y) i — > Fy(Tf o t =+oo). 
We start with a general formula which gives the density of the r.v. Gf. 
Theorem 1. 

Assume that H is of C 1,2 class on ]0, £(/)[x]£, +oo[ and is such that, for every x £]£, +oo[ and t > 0: 

llm a \a 9{t,x,y) = 0. 
y^+oo s'{y)oy 

Then, the density of the r.v. Gf under ¥ x is given by: 

x{ f ] ~ ~*UW) dy— ly=m 

Note that for x > f(0), this density is defective. 



dt 



(o<t< C(/))- 



(4) 



(5) 
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Proof. 

From Assumption (3), it is clear that Gf < ((f) a.s. Let < t < ((f): 
P* (G f < t) = F x (G f < t n X t > f(t)) 



+00 

fit) 

+OO 

fit) 

+OO 

fit) 



Px (Gf < t\X t = y) q(t, x, y)m(dy) 



F y (Tf a g t = +00) q(t, x, y)m(dy) (from the Markov property) 



H(t,y)q(t,x,y)m(dy). 



(6) 



Observe that, still from the Markov property: 

F y (Tf = +oo\F t ) = F Xt (Tf o9t = +00) l {Tf>t} = H(tAT f ,X tATf ). (7) 

Set (Z t = s(Xt /K T e ),t < ((f))- By construction of the scale function s, the process Z is a continuous local 
martingale. But, for X > /(0), T f A((f) < T e /\((f), hence the process (X tATf = s _1 (Z MT/ ),t < ((f)) is 
a semimartingale. Therefore, applying Ito's formula to (7), we deduce that the term with hnite variation 
vanishes, i.e. H is a solution to the partial differential equation: 



dH 

QH + — = on the domain {(t, y)\ y > f(t)}. 



Now, let us differentiate (6) with respect to t 

dH(t,y) 



? x (Gf edt)= I 
Jf(t 



at 



q(t,x,y)m(dy) 



+ °° tth ^g(t,x,y) 

H (^y) — ^7 — m ( d v) 
fit) m 



(8) 



(9) 



since H (t, f(t)) = 0. Integrating by part the second integral, we obtain: 



+ °° , dg(t,x,y) 

H (^y) qT rn(dy) = 

fit) dt 



H(t,y)Gq(t,x,y)m(dy) (from (1)) 

/(*) 

H (^y) 5 d V 

fit) 

+°° 8H(t,y) dq(t,x,y) 



H(t,y) 



dys'(y)dy 
J /(t) 

dH(t,y) dq(t,x,y) 



f(t) dy s'(y)dy 
dy 



dy 



hit) s '(y) d v d v 

(since H is bounded and +00 is a natural boundary, sec [BS02, p. 20]) 



9H(t,y) 



s'(y)dy q ( t,x,y ^ 

g(t,x,f(t))dH(t,y) 
s'(f(t)) dy 



d 2 H(t,y) 
fit) ' hit) 9ys'(y)dy 

+ OO 



q(t,x,y)dy 



\y=f(t) 



QH(t,y)q(t,x,y)m(dy). 



fit) 



It only remains to plug this relation in (9) and to use (8) to get the desired result. 
Example 2 (Transient diffusions and constant boundaries). 

Let (X t , t > 0) be a transient diffusion whose scale function is of C 2 class, with the normalization s(+oo) 
0. Let a > £ and choose f(t) = a. Then, for a < y < b, we have, by definition of the scale function s: 



□ 



(T b < T a ) 



s(y) - s(a) 
s(b) — s(a) 
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and, letting b — > +00, since +00 is a natural boundary: 

P„ (T„ = +00) = 1 - 44 
s(a) 

so we recover the well-known formula: 

P x (G a g dt) = - 9 ^f\ a ^ db, (0<t< +00) 
s(a) 

see Pitman & Yor [PY81]. Note that in this case, Equation (4) reduces to lim q(t, x, y) = 0, which is of 

y— > +00 

course always satisfied in our set-up. 

Remark 3. More generally, it can be proven that, for monotone functions /, Equation (4) is automatically 
satisfied, provided that H is smooth enough. 

• Assume first that / is increasing. Let t be fixed and set a = f(t). Then, for y > a, by the continuity of 
paths and the strong Markov property: 

T v {T a = +00) = P„(T /o9t = +00) + P y {T a = +c» n Tf g t < +00) 

r+oo 

= H(t, y) + / V y (T a = +(x\T fo6t = s)¥ y (T fo9t e da) 
Jo 
r+00 

= H(t,y) + / P/(t +s) (T a = +oo)F y {T foet £ da) 
Jo 

(see Peskir [Pes02]) 
= H(t,y) +E„[tt(T /oflt )] where *(s) = F f(t+s) (T a = +00). 

Observe now that, since 4" is increasing, all three functions are increasing functions of y. Therefore, 
we deduce, with the same normalization as in the previous example (/ is increasing hence X must be 
transient), that: 

Q < dH(t,y) < d¥ y (T a = +00) _ 1 



s'(y)dy s'(y)dy s(o)' 

which implies, since lim q(t,x,y) = 0, that Equation (4) is satisfied. 

y— > +00 

• Assume now that / is decreasing. In this case, the function t 1 — > H(t,y) is increasing, hence, from 
Equation (8): 

d 2 H{t lV ) _ dH(t,y) 
dys'(y)dy P[V> M ' 

we deduce that the function y 1 — > — ttt^ i s a positive and decreasing function. Therefore it is 

s'{y)oy 

bounded at +00 and Equation (4) is satisfied. 

We now give a few examples involving Brownian motion and Bcssel processes. From Theorem 7, we 
only need to compute P y (Tf g t = +00) to obtain the density of the last passage time to the boundary /. 

Example 4 (Bessel processes and straight lines: f(t) = a — bt). 

Let (Rt,t > 0) be a Bessel process of index v > — 1. (R t , t > 0) is a diffusion whose generator reads: 

gM = l g , 2„ + l 9 



2 dx 2 2a; 9x ' 
Its speed measure and scale function s^ v > are given by: 

fmW(dx) = 2o; 2 ' y+1 d2:, 
1 ( s M)'(x) =^ 2,y - 1 . 
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We denote by the law of (R t ,t>0) when started at x and by its transition density function with 
respect to nv- v ': 



with I v the modified Bcssel function of the third kind. 

Choose f(t) = a — bt with a, b > 0, and thus £(/) = — . Letting A tend toward in Theorem 5.1 of Alili 

a 

& Patie [AP10], we deduce that, for y > a: 

< + ~> - - (ry - £ r t^ M (s- •) 

where if,, denotes the McDonald function of index v. Then, since: 

®%\T Mt < +oo) = PM(T B _ W _ 6 . < +cx)) 

we obtain : 



8y lv=a - M 
and finally, for < t < — : 



.V^ i , \ 1, _„ / x 2 + (a-6t) 2 \ T (x[a-bt)\ dH(t,y). 
'>(<?._*. e A) = -(a - &t) + a: "cxp ^ L_ J /„ ^ j 



f/f. 



Example 5 (Reflected Brownian motion and straight lines: f(t) = a — bt). 

In particular, when v = — \ (i.e. when (X t ,t > 0) is a Brownian motion reflected at 0), then K_i = Ki 
and the above formula simplifies to: 



I y 27r(a - 6t) IV 2 \ ) Isfhrt 

2 r _ 2 

where Erf denotes the error function: Erf(z) = —= I e x dx. 

V* Jo 



Example 6 (Bcssel processes and squared root boundaries: /(f) = ay 1 + 2^ft ). 



Let (Rt,t > 0) be a Bessel process of index v > — 1 and choose f(t) = a\Jl + 2*yt with a > 0, 7 < 0, and 
thus C(/) = -— . Define 

r * = —z • 

2 7 

Then, it is known that the process (e _7 'i? T ( t ), t > 0) has the same law as a radial Ornstein-Uhlenbeck 
process with parameters v and 7, i.e. 

(X t = e-^R T(t) ,t>0) 

is a linear diffusion with characteristics: 

i^\dx) = 2a; 2 ' y+1 e^ :l:2 da;, 
r+00 

s^(x) = - / y- av - 1 e n dy. 
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Now, we obtain by time change, since r is increasing: 

inf (u > 0, R u = av /l + 27u| = r (inf {t > 0, R r{t) = ae 7 *}) 

= r(inf {t > 0, X t = a}) 

and we deduce from Example 2 that: 



p(")(T / = +oo) = l- V-fZT- 



Then, since 



(/ o 8 t ) (u) = a^l + 2 1 (t + u) = ay/1 + 2~/t Jl + ^^ u, 
Vi v) (Tf°6 t =+oo) = l 



we obtain: 

gK"' l+2-,t> (y\ 



2-1 

S 



Ktt^) (ay/1 + 2-yt) 



and finally, for < t < — -7- : 

2 7 



(G f G dt) = - - 2T r-g w (t,x,ay/l + ¥ft) dt. 



2.2 Upper boundaries 

We now study the case when the process (X t ,t > 0) remains under the boundary / after the last passage 
time Gf. Let £ £ [— oo,+oo[. We make the following assumption on the nature of the boundary point t 
(see [BS02, p.19-20]): 

• If £ = —00, we assume that £ is natural. 

• If £ > — 00, we assume that £ is entrance-not-exit. 
These hypotheses ensure that for every t > and x > £: 

dq(t,x,y) 

llm „ // \* = °- 
y->t s'(y)dy 

Let / : [0, +oo[ — >]t, +oo[ be a continuous function, which is of C 1 class on ]0, +00 [, and such that 
C(/) := inf{i > ; f(t) = £} = +00. We assume that : 



Vx el, ¥ x [ lim X t - f(t) < = 1. 

yt-y+oo J 

This implies that after Gf, the diffusion (X t ,t > 0) remains under the boundary /. 
Let H be the function defined by: 

H ■ (t, y) 1 — ► Fy(Tf o t =+00). 

Theorem 7. 

Assume that the function H is of C 1 ' 2 class on ]0, +oo[x]£, +oo[ and is such that: 

Kmj¥^q(t,x,y)=0. (10) 
v^e s'(y)dy 

Then, the density of the r.v. Gf under ¥ x is given by: 

Pl(G/efl) = _i2^^| >=/(()(it ((>0 >. (n) 
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Proof. 

The proof is of course very similar to the previous one. For < t < +00, we have: 

/■/(«) 

Vx(G f <t)= H(t,y)q(t,x,y)m(dy) 
Je 

with the function H solution of 

da 

QH+ — =0 on the domain {(t, y); y < f(t)}. (12) 
at 

(Note that, by hypothesis, Ti = +00 a.s. hence, since s is a function of C 2 class, X is a semimartingale) . 
Then, we make the same integrations by parts, noticing that the terms between brackets cancel since: 

doit x v) 

1) H is bounded and > 0. 

s'(y)dy y^t ' 

ii) lim — ^ '' ^ q(t, x, y) = by hypothesis. 
y^e s'(y)dy 



□ 



Example 8 (Bessel processes and straight lines: f(t) = a + bt). 

Let (Rt,t > 0) be a Bessel process of index v > and choose fit) = a + bt with a, b > 0. Then, from 
Theorem 5 of Alili & Patie [AP10], for y < a: 

n ) (T a+b .€du) = exp[-(y a ) -u 1 + -u) ^ exp - — — — \du, 



^ a 2 v J v +x{ju,k) \ 2a(a + bu)^ 



where J v denotes the Bessel function of the first kind and (j^k)k>o the ordered sequence of its positive 
zeroes. Then, 

H(t, y) = l- ¥^(T a+bt+b . < +00) 

and, for x > 0: 

m Ga+ , s *) = -i(a +6 *r^«p (- - 2+( ; +bt)2 ) /„ ^i^*, 

i(a + exp r- " 2 + ( ? + bt) i J* ( (-6 + ^(a + W)) * 



2V ' 1 V 2 * J \ t 

with 

^ ) = /o eXP ^T M Jl 1 + ^ U J g (« + bt)^J v+1 ^ k ) 6XP I," 2(a + bt)(a + b { t + u)) ) ^ 
Example 9 (Brownian motion and square boundary: f(t) = a + bt 2 ). 

Let (B tl t > 0) be a Brownian motion and define f(t) = a + bt 2 with a, b > 0. From Salminen [Sal88], we 
know that: 

My < a, Fy(l / e du) = 2(bc) exp « J ^ exp ^ — u j 1 



fc=0 



Ai'(A fc ) 



where Ai denotes the Airy function, (Xk)k>o its zeroes on the negative half-line and c = (2b 2 ) 3 . Now fix 
t > 0. Applying the Cameron-Martin formula, we deduce that: 

^ { y-uh T f e d«) = exp ^-2fa(a - y + ^ 2 ) - (2&t)^ 2 - Pj/-6t 2 ( T / € d«), 
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where denotes the law of a Brownian motion with drift fi started at x. But, since 

inf{u > 0; (y - hi 2 ) + B u - 2btu = a + bu 2 } = m£{u > 0; y + B u = a + b(u + t) 2 }, 

we obtain that F y (T fo6i < +00) = P ( ~_^(Tf < +00) and 



H(t, y) = 1 - 2{bc) 2 cxp ( ~b 2 t 3 - 2bt(a - y) 



r+ °° ( 2, 2/ , a \ ^ / Afe 
exp I --V(u + ty I 2^ ex P ( — ' 



I] 



Afe ^ Ai(A fc + 2bc(a - y + bt 2 )) 
Ai'(Afe) 



du. 



Finally, for 



\{G f G dt) 



1 



2V27rt 
1 

2\/2~7rt 



cxp 



(x-a-bt 2 ) 2 \ dH(t,y) 
2t ) dy 



\y=a+bt 2 dt, 



~ { - a 2t ht2r ){-2bt + ^{a + bi 2 



with 



+00 



My) = -2(6c) 2 cxp ( |6 2 t 3 ) / cxp ( --b 2 {u + tf ^exp ^ 



k=0 



Afe \ Ai(Afc + 2bc{a - y + bt 2 )) 



Ai'(Afe) 



du. 



3 Martingales methods 

We now present a method to obtain explicit expressions for the function H (associated with an a priori 
implicit boundary). Let £ <G [— oo,+oo[ and / : [0, +00 [1 — >]£, +oo[ be a continuous function. In this 
section, we shall restrict our attention to lower boundaries, and make the following Assumption: 

Assumption 10. If Q{ f) = +°°,- we assume that: 

lim X t = +00 a.s and P x [ lim X t - f(t) > ] = 1. 

t— »-+oo \ t— ^+00 J 

Consider the domains: 

2> = {(t,y)e]0,C(/)[x]4+oo[; y>f(t)} and dV := {(*, y) G]0, C(/)[xR +oo[; y = /(*)}. 

Lemma 11. Assume that Assumption 10 holds and that there exists a function H : T> 1 — > [0,1] of C 1,2 
class on T> which is solution of the following problem: 



OH 

QH + -7^- = 0, with the boundary condition H(t, y) = 1 on dT>, 



H is decreasing in y 

l/C(/)_=+°°, 

lim H(t, y)=0 

K t,y— j-+cxd 



(13) 



7/C(/) <+oo, 
Vy e]^,+oo[, lim = 0. 



Then, 



H{t,y)=F y (T Mt < +00). 
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Proof. 

Let t < C(/) be fixed. Applying Ito's formula, we deduce that the process (H(t + s, X s ), s < (£(/) — t) A Tf g t ] 
is a positive and continuous local martingale. Let s < £(/) — t be fixed. From Doob' stopping theorem, 
with f{t) < y < a: 



H(t, y) = E v [H(t + s A T a A T fo9t , X sATaAT/o6t ] 



EL 



H(t + T a , a)l {Ta<sAT/o9t} + P y (T /oflt < s A T a ) 



+ s,X s )l {s<ToAT/oe j 
We first let a — > +cx) to obtain, applying the dominated convergence theorem: 

H(t,y)=E a [H(t + S ,X s )l {s<Tfoet} ] +¥ y (T fo g t < s) 
Then, we must distinguish between two cases. 

1. If £(/) = +00, then, we assumed that lim X s = +00 a.s., so: 

— > a.s. 



H(t + s,X s ) - 

s— ^+00 

and the result follows from the dominated convergence theorem. 
2. If C(/) < +00, then, since F y (X c(f) _ t = I) = 0: 

lim H(t + s,X s ) = a.s. 
»^C(/)-i 

and the result follows once again from the dominated convergence theorem. 



We shall now give some examples of functions H which are solutions of this problem. 



□ 



3.1 Martingales constructed on the resolvant 

We follow the idea of Robbins & Siegmund [RS70]. Let us define, for A > 0, the resolvent kernel of 
{X t ,t> 0) by (see [BS02, p.19]): 

u\(x,y) = / e~ xt q(t,x,y)dt. 
Jo 

Let F be a finite measure on [0, +00 [, and define: 

r+00 

H(t,y)= / e- xt u x (0,y)F(dX). 
Jo 

We assume furthermore that: 

1. if I = — 00, then £ is a natural boundary, 

2. if I > —00, then t is an entrance- not-exit boundary. 

This implies that limuA(0,y) = +00. Now, since the function y 1 — > u\(0,y) is strictly decreasing, we 
may define a function / : [0, +oo[ — >]£, +oo[ by 



Vt>0, H(t,f(t)) = l, 



which is such that £(/) = +00. Then: 
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Proposition 12. The density of Gf is given by: 

Proof. We need to check that the function H satisfies the hypotheses of Lemma 11. 

By construction, H is a solution to the partial differential equation (13). It is also a decreasing function of y 
and, since +oo is a natural boundary, we have from the monotone convergence theorem, lim H(t, y) = 

t,y— >+oo 

0. 

□ 

Example 13 (Brownian motion with drift). 

Let (Xt,t > 0) be a Brownian motion with drift (i, and choose F(dX) = 2fiaSo(dX) + 2y26 + /j, 2 <5&(<fA). 
Then: 

H(t,y) = ae- 2 ^ + ^-(V^+^y _ 

Define: 

( ^2b + n 2 + p) y + ln( 1 - ae~ 2 ^) 
V{y) = ~- — b and f(t) = ip- 1 (t). 

f is a decreasing function such that lim f(t) = — Then: 

t— S-+oo 2\X 



dt. 



Example 14 (Bcssel process of dimension 3). 

Let (X t ,t > 0) be a Bessel process of dimension 3, and choose F(dX) = a5o(dX) + 8b(dX). Then: 

-bt-V2by 
H(t,y) = * + P - . 

y y 

Define: 



V2by + ln(y - a) _ 1 
(p(y) = and f{t) = <p (t). 

f is a decreasing function such that lim f(t) = a. Then : 

3.2 Martingales constructed on the transition density 

We assume in this subsection that (X t ,t > 0) is defined on (0, +oo[. 

3.2.1 Preliminaries 

We start by recalling a few properties of the transition density (t,x) i — > q(t,x,Q). From Kotani & 
Watanabc [KW82], it is known that : 



y p(z)s'(z) \ _ 



ii 



In [KW82], this formula was obtained in the special case s(x) = x, but, assuming that s is a strictly 
increasing function of C 1 class, Formula (14) follows easily from the fact that s(X) is a diffusion on 
natural scale. In particular, for x > 0, we deduce that: 

lim q(t, x, 0) = 0. (15) 

t-»o 

Next, we define a new diffusion (X t , t > 0) whose speed measure m(dx) = ~p(x)dx and scale function s are 
given by Bianc's transform: 

'-p{x) = (m{[0,x])) 2 s'{x) 

- (x) = __L L_ 

S[X) m([0,+oo[) m{[0,x])' 

It is known from [PRY10, Chapter 8] that the transition densities of X and X satisfy the following relation: 

/>oo 

q(t,y,0)= q(t,0,z)m([0,z])s'(z)dz, (16) 



which implies in particular that for every t > 0, the function y i — > q(t, y, 0) is decreasing and tends toward 
as y — > +oo. 

3.2.2 First example 

Let £, c > and consider the function: 

H(t,y) = -q(C-t,y,0). 
c 

When c < inf q(( — t, 0, 0) we may define a boundary / by: 

q(C-t,f(t),0) = c. 
In this set-up, C = inf{t > 0; f(t) = 0} = C(/) from (15). Then: 
Proposition 15. The density of Gf is given by: 

P x {G f edt)- — cjfffty ~d y \v=fW dt [0<t<Q. 

Proof. From (1), the function H is a solution of the PDE (13), which is decreasing from (16) and, for 
every y > 0, lim H(t, y) = from (15). 

□ 

Note that the martingale (H(t, X t ),t < £) appears as a density when constructing diffusion bridges 
via Doob's /i-transform, see Fitzsimmons, Pitman & Yor [FPY93]. 

Example 16 (Radial Ornstcin-Uhlcnbeck processes). 

Let (X t ,t > 0) be a squared radial Ornstein-Uhlenbeck process with parameters v > — 1 and 7 < 0, see 
Example 6. Its transition density function reads (see [BS02, p. 142]): 

Q V, V, ^ 2 ^T(v + 1) (sinh(7i))^ +1 P ^ 2 sinh( 7 i) J ' 
Let a €]0, 1[, C S]0, +00 [ and choose: 

C = a^(C,0,0) 
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With these values, the boundary /(">t) is then defined, for t < £ by: 



2 s _ u , , A f^Mi(C-t))^ 



sinh( 7 (C - t))eT(C-*) In U / 7 ^ x ;; e 7(^+i)t 

\ 7 \ V sinh(7C) / y 

and the density of the r.v. G/ is given by: 

P,(G /( „ 7) e di) = P^t^t^ (/^(t)) 2l ' +2 e-T( /< '" 7) W) 2 9 ^(t,a :) /^)(t))^ (0 < i < C). 
J csmh(7(4 — t))\ / 

Example 17 (Bessel process). 

Letting 7 — > in the previous example, we obtain the following boundary for the Bessel process of index 
v. 



^ 2(C-t)ln^(^l-- 



and the density for the last passage time to f^ u ' reads: 

PM(G /M e dt) = (/W^J^^M^xjMW)*, (0 < t < 0, 

• , a 

With C = — ; r — . 

3.2.3 Second example 

We assume that is an entrance- not-exit boundary point and that (X t , t > 0) is transient and goes towards 
+00 as t — >• +00. Let h : [0,+oo[ — >-]0,+oo[ be a bounded function of C 1 class with bounded derivative. 
We consider the function: 

r+00 

H(t,y)= h(t + u)q(u,y,0)du. 

Jo 

Since we have assumed that is entrance-not-exit, U\(0, y) > +00, i.e. q(t, 0, 0) is not intcgrable at 0. 

Therefore, for every t > 0, the function y 1 — >• H(t,y) is decreasing from +cx) to 0, and we may define a 
function / by: 

f +00 

h(t + u)q(uj(t),0)du = 1. 

Note that by construction, £(/) = +00. 
Proposition 18. TTie density of Gf is given by: 

'AO, e it) - -Sfci^a (/ + ~ „(„ + fl^U,,,,*) « > 0). 

Remark 19. We may recover (partly) the result of Subsection 3.1 by taking h(u) = J +O ° e~ Xu F(d\). 
Indeed, in this case, from Fubini-Tonelli: 

+00 r-\-oo / />+oo \ />+oo 

h(u + t)q(u,y,0)du = J (j e- x{t+u) F(d\)\ q(u, y, 0)du = J e~ xt u x {y,0)F(d\). 
Proof. We need to prove that H satisfies the hypotheses of Lemma 11. 

Since h is bounded, and for y > 0, u 1 — > q(u, y, 0) is integrable, the dominated convergence theorem 
implies that : 

lim H(t, y) = 0. 

t,y— >+oo 
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Then, integrating by parts, for y > 0: 



dH(t,y) 
dt 



r+oo 

/ h'(t + u)q(u, y, 0)du (which is finite since h! is bounded) 
Jo 



+00 



+00 



dq(u 7 y,0) 



h(t + u)q(u,y,0) — / h(t + u) ^ — - — du 

Jo Jo du 

+00 

h(t + u)Qq(u 1 y, 0)du 



-QH(t,y), 



which ends the proof. 



□ 



Remark 20. We may remove the hypothesis (X t , t > 0) is transient if we replace the assumption on h by : 
h is a decreasing and intcgrable function of C 1 class. Indeed, since for y > the function u 1 — > q(u, y, 0) is 
bounded, H is well-defined and so is L 00 h'(t + u)q(u, y, 0)du. Besides, we still have lim H(t, y) = 

from monotone convergence. 

4 Inverting time 

Consider a diffusion (X t , t > 0) enjoying the inversion property in the sense of Watanabc [Wat75], i.e. such 
that the process (Xt = tXi , t > 0) is also a linear regular conservative diffusion. Let / be a continuous 
function and define: 

7 = *•—>*/ (T 



Then, if X ^ lim/(t) 



t->o 



inf {t > 0; X t = /(*)} = inf It > 0; tXi = tf Q 



inf U > 0; Ii = / 



sup{u > 0; X„ = /(«)} ' 

In particular, under the hypotheses of Theorem 1 or 7, the density of the first hitting time of / by 
(X t ,t >0) admits the expression: 



dt. 



4.1 Brownian motion 



Consider a Brownian motion (B t ,t > 0) started from 0. Then, it is well-known that (B t ,t > 0) is also a 
Brownian motion started at 0. 

Example 21. Take f(u) = a + bu 2 with a,b> 0. Then, J(t) = at+j and : 

inf it > 0; B~t = at 



t J sup{u > 0; B u = a + bu 2 } 
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1\ 1 If 6 X 2 



From Example 9, we deduce that: 

¥oiT f edt) = *{l)^ cxp {-Yt{ at + V ]dt 

where is defined by : 



with 



/ 2 



^( y ) = -2(6c) 2 exp(|U / exp £ e xp(^ 



3*V io " V 3 \ tj ) £j " V c 7 Ai'(A fe ) 



with (Afc)fe>o the negative zeroes of the Airy function Ai and c = (2b 2 ) 1 ' 3 . 
4.2 Bessel processes with drift 

Consider a Bessel process (Rt,t > 0) with index > and drift c > started from a:. (Rt,t > 0) is a 
diffusion whose generator is given by: 



2 dx 2 \ 2x I v (cx) J dx 

This process was first introduced by Watanabe in [Wat75] as a generalization of Bessel processes. For 
integer dimension of n = 2(y + 1) e N, (Rt,t > 0) has the same law as the norm \\B t + /?. t]| where 
(B t ,t > 0) is an n-dimensional Brownian motion and = c. 
The scale function of i? is given by: 

1 MJ (r(^ + l)) 2 V2J xl 2 {cx) 

and its speed measure by: 

2i' 



mS v > c \dx) = 2(T(v + l)) 2 Q xl 2 (cx)dx 



From Watanabe [Wat 75, Theorem 2.1], the process \^R t — tRi,t > 0J is a Bessel process with index v 
and drift a; started from c. 

Example 22. Take f(u) = a + bu with a, 6 > 0. Then f(t) = b + at and 

1 



inf {t > 0; i? t = b + at} 



sup{u > 0; R u = a + bu} 
which leads to the density of the r.v. TV|_ a . under the form: 

(T b+a . = (j) (* (a + ~) ) ~" exp (- ^ 2+ ^ + b)2 ) I„ (s(<rf + 6)) <K 
where the function $ is given by (see Example 8): 



•G)-M)~VH 

with 



ii 



where denotes the Bessel function of the first kind and (j v ,k)k>a the ordered sequence of its positive 
zeroes. 
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5 On an integral equation 

We set: 

Integrating (5) and (11) with respect to t, we deduce the following corollary: 

Corollary 23. Assume that the hypotheses of Theorem 1, resp. 7, are satisfied. Then, the function 4> is 
a solution of the following Fredholm equation of the first kind: 



• For lower boundaries: 

K(/) 



$(t)q(t,x,f(t))dt = l, x<f(0). 



Resp, for upper boundaries: 

r+oo 

$(t)q(t,x,f(t))dt = l, x>f(0). 

In some particular cases, this equation may charaterize uniquely $, hence the law of Gf. We give 
below a few examples of this situation, where time inversion is involved. 

5.1 A link with time inversion 
Theorem 24. 

i) Let (Bf\t > 0) be a Brownian motion with drift ji > started from x and f be a continuous function 



[0, +oo[ such that 



Assume that the equation 



/(0) = and lim ^ < fj,. 

t— >+oo t 



/■+00 

Vx<0, / $(t)q(t,x,f(t))dt = l, 
Jo 

admits a unique solution $. Then : 

Pj>° (Gf G dt) = <S>(t)q(t,x,f(t))dt. 

ii) Let (Rt,t > 0) be a Bessel process with index v > —1 and drift c > started from x > 0, and let f 
be a continuous function on [0, +oo[ such that: 

/(0) = and lim ^ > c. 

t— >+oo t 



Assume that the equation 



r+oo 

Va->0, / m)q^ c \t,x,f(t))dt=l, 



admits a unique solution $. Then for x > 0: 

P^ c) (Gf G dt) = <$>(t)q^- c \t,x,f{t))dt. 
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Proof. 

i) The transition density of (B^ ,t > 0) reads: 

1 



q(t,z,y) 



2V2TTt 



exp -fj,(z + y)--z-t- 



•it 



Observe first that: 

q(t, z, y) = q(t, 0, y) exp ( -— + y - A 
hence, by hypothesis, $ is the unique solution of the equation: 

' z 2 , zf(t) 



Vz < 



0, J $(t)q(t,0,f(t))exvl-- 



or, with the change of variable s = j: 



Vz < 0, 



+oo / 2 



( -,o,/ 



exp 



t 



■s + zsf 



dt = e" z , 



1\\ ds 



a fJ.z 



Now, let /(,s) = sf ( — ) and consider, for z < 0, the exponential Brownian martingale ^cxp ^z_B s — -ys j , s 



Doob' stopping theorem implies, since fi > lim /(s) by hypothesis: 

s— >0 



exp ( zB T j - —T- f 



+00 



(T 7 e ds) 



Since / is a positive function, Tj admits a density and by unicity of the solution $, we have: 



Next, the Cameron-Martin formula gives: 

(r 7 < *) = e 

i.e. the density of the r.v. Tj equals: 



-,o,/ 



<is. 



exp I nBr 7 - y Ij - I 1{t-<s} 



f> (Tj e ds) = exp - ^s - fix) 1$ ( -) q 



1 



: ,0,/ 



v 

c- 



1 



ds. 



The result then follows from the time inversion property. 

it) We shall proceed similarly for Bessel processes with drift. The transition density of (Rt,t > 0) reads: 

i fc\*» cx p(^4^) 



(t,z,y) 



2t(r(i/+l)) 2 V27 I u (cz)I v {cy) 
hence, let $ be the unique solution of the equation: 



exp 



z 2 + y 2 



2/ 



(f) 



Vz > 0, 



+00 



*(*) 



o 2i(r^ + l))2 V2 



xp / c 2 t (/(t)) 2 \ 



Iu{cf{t)) 



2\ 



ex P — 777 



2i / I v {cz) 



Idt = 1. 
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Consider the local martingale under 



'm^.) ■' ao J- 

Applying Doob' stopping theorem to (MT-AT e At,t > 0) and letting t — > +00 and e — >• 0, we obtain (since 
c < lim f (t) and /(£) = o(t) by hypothesis): 

t— >0 +00 



I v (xc) 

Now, from Watanabc [Wat 75]: 



\ e~ xt yJ — _ >-^' x) (T T £ dt) . 

Jo Iu{xf(t)) ' f 



{",*) ( lim *L = x \ = l, 

t-S.+oo £ 



hence, since /(f) = o(t), we deduce that, for x > 0, the r.v. Xy admits a density (see also [Leh02]) and, 
setting 2A + x 2 = z 2 , we obtain from the unicity of the solution $: 

WW) Fc l J / ed V-^ 2(r( , + 1))2 UJ I,(c/(i)) 
The result then follows once again from a time inversion argument. 

Observe that, if Pc (Ty < +00) = 1, (for instance if / is of C 1 class in the neighborhood of 0), then the 
result also holds for x — 0. 

□ 

It might be noticed that for these two processes, the proof above shows anew the phenomenon of separation 
of variables which appears in the law of Gf. We shall now give an example of each situation. 

5.2 Brownian motion with drift and f(t) = a + b\fi 

Let (Bt\t > 0) be a Brownian motion with drift /1 > and choose f{t) = a + b\/t with a, b £ R. 
Proposition 25. The density of the last passage time G a+b ^ is given, for every x £ R, by: 

(G a+b ^. € dt) = cxp (-^^ + - M* -»)-?) * 

where the function if has Mellin's transform: 

+00 ^ 
t x ~ 1 Lp(t)dt = 







In particular, if x = a, we deduce that: 

■(-*) 



exp 

Ea ' 



M^-^-aA+i (6)' 



Proof. We need to prove that the equation: 



/■+°° ( z 2 zb\ 

Vz > 0, y $(t)g(i, 0, exp [- — - — J = cxp(-^z + 2^a) , (17) 
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admits a unique solution $, in order to apply Theorem 24. Let us recall the following formula ([GR07, 
3.462-1]): 

z v ~ x exp [ — ^-z 2 — jz) dz = [ L exp [ -73 J D_ v ( — 7= ] with f3, v > 0, and 7 £ R, 



where D- v denotes the parabolic cylinder function of index —v. Integrating Equation (17) with respect 
to z u ~ 1 dz and applying Fubini-Tonelli, we obtain: 

r+00 /h 2 \ r(v) 

/ $>(t)q(t, 0, bVt)Y(v)t v/2 exp — D_ v (b) dt = exp (2/m) 
Jo V 4 / 

which, by setting A — 1 = ^, gives the Mellin's transform: 

+00 exp ( 2\ia — ^j- 

<5>{t)q(t,0,bVi)t x - 1 dt 



A ^- 1 ^_ 2A+1 (6)' 

Since Mellin's transform is injective, the result follows from Theorem 24. 



□ 



5.3 Bessel process with drift and f(t) = \J at 2 + bt 

Let (Rt,t > 0) be a Bessel process with index v > — 1 and drift c > started from x, and choose 
/(t) = \Jat 2 + bt with b > and -^a > c. 

Proposition 26. XTie density of the last passage time G ^/ a 2 +b- is given, for every x > 0, by: 

* V Va - +6 - 7 V V at) ) ct^aWVbt L(cx) "y t J 
where the function if has Laplace transform: 

+» At _ exp(^)M_ AW2 (^) 
exp (|) M_ Ail//2 (2) 
Proof. We need to prove that the equation : 

$(*) e -4* / x 2 + (/(t)) 2 ^ / W^+to ^ 



Vz > 0, 



o 2t(r(i/+l))» V2V J„(c/(t)) 

(19) 

admits a unique solution $ in order to apply Theorem 24. To simplify the notation, set 

$(f) (c^ e"# / (f(t)) 2 



* {t) = 2t(r(v + i))2 UJ /Jcf(O) cxp 



2t(I> + l)) 2 V2^ J„(c/(t)) 2i 
Recall the formula [GR07, 6.643-2], for a > 0,f3 e R and A + f + \ > 0: 



where A/_a.„ denotes the Whittakcr function. We integrate (19) with respect to x e 2bX dx, to obtain 
2t ei (b\ , 2 f2b\ X /c 2 6\ ( c 2 b\ 



'O ga^ -V/. ( 5 J * = - ("J exp ^- j M_ A ,„ /2 v 2o 
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This expression simplifies to: 



n 



t fc 2 b b\ m -X,u/2 I Ti£ 



+00 I 7 -1 -1 / J. >\ M. 



at + 6 (1 + A)A c V M_ A>I//2 (§) 



Finally, we make the change of variable e u = 1 H : 

at 



^ AU * ( " 7\ ) „, du = ^ exp 



o \a(e"-l)J (e«-l)2 6c ^ V 4a V M_ A ,„ /a $ 

which gives the Laplace transform of VP up to a few transformations. 

Remark 27. Taking A = a — i ^r L , we obtain, since y/a > c: 



□ 



I 



+00 eX p^W ,„+i /9 f#l 



exp(|)M_ Q+ ^ W2 (|) V c ^ VP J 



where Qi 1/ ' 7 ' ) denotes the law of a radial Ornstein-Uhlenbeck process with parameters v and 7 started at 
x. Therefore: 

<PW = - P f-^±±t) (£) ^ (T^ e dt) . 



Corollary 28. Ta&e f(u) = ^fwu? + bu with b > and -^/a > c. TTien /(t) = \/a + bt and: 
inf{t >0;R t = Va + bt} 



sup{w > 0; R u ~ a/ an 2 + bu} 
which leads to the density of the r.v. T v / a+b- under the form: 

Pi^PV^r G = V fin fl + - 7 1 , ^/r^/, fi-^+w) dt 

v + V V aJJ cy^Tbi Iv{cx) V / 

where the function ip is given by (18). 
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